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, |_H ' Abstract. In this paper we investigate the Sunyaev-Zel'dovich (SZ) effect and the X-ray surface 



t 



brightness for clusters of galaxies with a non-spherical mass distribution. In particular, we con- 
sider the influence of the shape and the finite extension of a cluster, as well as of a polytropic 
thermal profile on the Compton parameter, the X-ray surface brightness and on the determina- 
tion of the Hubble constant. We find that the non-inclusion of such effects can induce errors up 
to 30% in the various parameters and in particular on the Hubble constant value, when compared 
with results obtained under the isothermal, infinitely extended and spherical shape assumptions. 
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1. Introduction 

Over the last few years, studies on galaxy clusters using X-ray emission observations have been 
a source of a tremendous increase in the litterature, especially those using Sunyaev-Zel'dovich 
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1995, Birkinshaw 



1999) is one 



(SZ) effect. The SZ effect (Sunyaev & Zerdovich|197g Rephaeli 

of the major sources of secondary anisotropics of the Cosmic Microwave Background (CMB) 

arising from inverse Compton scattering of the microwave photons by hot electrons in clusters 

of galaxies. 

Many different works have been developed during recent years leading to the use of this effect 



for studies of cosmology (Bernstein & Dodelson 1990; Aghanim, de Luca, Bouchet et al. 1997; 
Barbosa, Bartlett, Blanchard et al. 



1996 



Bartlett, Blanchard & Barbosa 



1998; Cooray 1998). 



Observations in the millimetre and submillimetre wavebands (Perrenod & Lada 1979; Chase et 



al. |1987t Silverberg et al. |1997[ ) give important information on the characteristics of clusters of 
galaxies. For example, by combining the SZ intensity change and the X-ray emission obser- 
vations, and solving for the number density distribution of electrons responsible for both these 
effects (after assuming a certain geometrical shape), the angular diameter distance to galaxy 
clusters can be derived. Assuming a cosmological model, this leads to an estimate of the Hubble 



constant Hq (Holzapfel, Arnaud, Ade et al., |1997| for Abell 2163; Birkinshaw & Hughes, |1994 
for Abell 2218). 

The SZ effect thus offers the possibility to put important constraints on the cosmological models. 
For this reason, different projects to measure the SZ effect are under way for example the MITO 



instrument (De Petris, Aquilina, Canonico et al., 1996) or (the longer term) the Planck mission 



(ESA report |1997| ). 

The SZ effect is difficult to measure accurately, since systematic eiTors can be significant. For 



instance, Inagaki, Suginohara & Suto (1995 1 made an analysis of the reliability of the Hubble 
constant determination based on the SZ effect. 

An additional effect arises if the cluster has a peculiar velocity (kinematic effect). Several papers 
discussed the influence of the kinematic effect on the measurement of the thermal SZ effect (De 



Luca, Desert & Puget, 1995 ; Audit & Simmons, 1999 for transverse clusters velocities). Note 
that the kinematic effect can thus be used to infer the peculiar velocity of clusters of galaxies. 



if the value of the Hubble constant is known (Rephaeli & Lahav, 1991 ; Haehnelt & Tegmark, 



1996). Another possible distortion on the SZ effect is due to gravitational lensing (Blain, 1998; 



Roettiger et al., |1997[ ). 

The extension and the geometry of hot gas distribution in clusters of galaxies is also an impor- 



tant source of systematic errors in the SZ effect. Cooray (1998 ) showed that projection effects of 
clusters can affect the calculations of the Hubble constant and the gas mass fraction. Recently, 



SuUcanen (1999) showed that galaxy cluster shapes can produce systematic errors in Hq mea- 
sured via the SZ effect. It is thus necessary to know, at least approximately, the shape of the 



clusters, for instance, if they are oblate or prolate (Cooray 2000 Hughes & Birkinshaw 1 998 ) or 
have a more general geometry. 



The /3-model (Cavahere & Fusco-Femiano 1976 ) is widely used in X-ray astronomy to parametrize 
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the gas density profile in clusters of galaxies by fitting their surface brightness profile. Neverthe- 
less, fitting an aspherical distribution with a spherical /3-model can lead to an important inaccu- 



racy (see Inagaki, Suginohara & Suto 1995 ). 

The aim of this paper is to investigate the influence of the shape and the finite extension of an 
ellipsoidal cluster gas distribution on the SZ effect, and to discuss the possible errors induced in 
the inferred value for Hq. The paper is organized as follows: 

In Section we present the calculations of the SZ effect and the X-ray surface brightness for an 
ellipsoidal shape with an isothermal profile and a finite cluster extension. Details of the calcula- 
tions are reported in two appendices. 

Section || is then devoted to a quantitative discussion of the incidence of these effects on the 
SZ measurements, in particular, the finite extension and the geometry (prolate and oblate) of the 
cluster. The influence of a polytropic thermal profile on the SZ measurements is also considered. 
The discussion and conclusion are given in Section Q. 



2. Basic equations of the SZ effect and the X-ray surface brightness for an ellipsoidal 
geometry 

Different X-ray surface brightness measurements in clusters of galaxies clearly indicate an as- 



phericity of the cluster shape. Fabricant, Rybicki & Gorenstein (1984) showed a pronounced 
ellipticity for the cluster Abell 2256 which indicates that the underlying density profile is aspher- 



ical. Allen et al. (1993) reached the same conclusion for the profile of Abell 478, and Hughes, 



Gorenstein & Fabricant ( 1988 ) for the Coma cluster. It is thus of relevance to study the influence 
of non-spherical shape on the results of clusters which have been reported so far. 



Given the above results, we have assumed an ellipsoidal /3 model: 

'3/3/2 



'^eyx-) ^y-] '^z) — ^e 



,p2 rj.^ 2 

i + li + ^ + i^ 



(1) 



Cl '^2 S3 

where rieo is the electron number density at the center of the cluster and /3 is a free fitting param- 
eter, which lies in the range 1/2 < /3 < 1- The set of coordinates r^, Vy and r^, as well as the half 
axes of the ellipsoid, (^i, C2 and (^3, are defined in units of the core radius re of the corresponding 
spherical shape. 

The fractional temperature decrement ATsz of the cosmic microwave background due to the 
SZ effect is expressed as 

(2) 



ATsz = 
with 



-/My 



1 



- 1 



where lu 



hv 



ksTcMB 



and 



(3) 



Tqmb is the temperature of the cosmic background radiation at z ~ (Tcmb ~ 2.728±0.002 



K, Fixsen et al. 1996) and ks the Boltzmann constant. The Compton parameter y is defined as 



2rr 



WeC2 



-arriedry 



(4) 
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We have chosen the line of sight to be along the r.y axis. I is the maximal extension of the hot 
gas in units of the core radius Tc along the line of sight, Te the temperature of the hot gas, me the 
electron mass, ctt the Thomson cross section and c the speed of light. 
The X-ray surface brightness of a cluster is given by: 

where z is the redshift of the cluster, which takes into account the cosmological transformation of 
the spectral surface brightness, ex the spectral emissivity of the gas, which can be approximated 



by a typical value (for Te > 3 x 10^ K, Sarazin |1986| ): 

ex = e Vr'e with e « 3.0 x 10"^^ n^ erg cm"^ s"^ , (6) 

where rip denotes the proton number density. The thermal SZ effect and the X-ray surface bright- 
ness depend on the temperature profile and of course on the density profile. We consider in the 
following an isothermal profile. The influence of a polytropic thermal profile is discussed in 
Section 3. 



2.1. Isothermal profile 



In this case (with Te = Teo) the Compton parameter and the surface brightness depend only on 



the density profile 

kBTeoUeorcCrT 



TTleC^ 



Isz with Isz = 2 / —^dr. 
ri, 



S 



X 



■Is^ With Is^ 



y ■> 

2 



dr„ 



(7) 



(8) 



47r(l + zy- JO 

where we introduce the structure integrals Isz and Isx ^ which depend only on the geometry and 

the extension / of the cluster along the line of sight Vy . 

With the structure integrals calculated in Appendix A, we obtain for y and Sx '■ 

, kBT,oaTneoC2r. ( rg ^^ '^^^^ 



Sxirx,rz) 



V2 2'2 

enl^y/TZC2rc 
47r(l + z)3 



B„ 



2' 2'2 



Cf 



^3 



-3/3+i 



B[W-l,^-B„Ail3- ^ 



2' 2/ ™ y f^ 2' 2^ 

where the cut-off parameter m, which depends also on the extension I of the cluster, is: 

777 ^ 

For the ratio between y"^ and S'x(?'i;, ?'z) we thus obtain 



(9) 

(10) 
(11) 
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\Ba 
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(12) 
(13) 
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3. Errors obtained in the quantities y, Sx and Hq 

3.1. Finite extension of clusters 

Since the hot gas in a real cluster has a finite extension, each of the observed quantities of the 
Compton parameter and X-ray surface brightness will be smaller than that estimated based on 
the assumption that I -^ cx). The incomplete beta-function Bm in Eqs. (y) and (jl^) can be seen 
as a correction term due to the finite extension, which together with the geometry of the cluster, 
enters through the m-parameter (see Appendices). As an illustration we report here the analysis 
of the influence of this correction for the simplest cluster case: isothermal (3 — 2/3-model with 
a spherical density profile (i.e. Ci = C2 = Cs ^ 1)' ^"d a line of sight going through the cluster 
center (i.e. Vx = r^ ^ 0). The reason for this choice is to be able to neglect, in this case, 
geometrical effects and thus to focus only on the modification due to the finite extension. 



0.20 



0.15 - 



0.10 - 



0.05 



0.00 




4 6 

Cluster extension 1 (in units of r^) 

Fig. 1. Relation between the ?7i-parameter of the incomplete beta-function Bm and the cluster 
extension I. The line of sight goes through the center of a spherical cluster (isothermal /3 = 2/3- 
model). 



In Fig. (|l|) we have plotted the parameter ttt, as a function of / for an isothermal (3 = 2/3- 
model. We notice that high values of I correspond to small values of m, which is almost vanishing 
for I > 6, and, therefore, the correction for the finite extension becomes negligible (to -^ for 
I -^ 00). 

We denote y(oo) and Sx{oo), respectively, as the expressions of the Compton parameter 
and of the X-ray surface brightness for a cluster with infinite extension, and y{l) and Sx{l) as 
the same expressions for a cluster with a finite extension / along the line of sight. The relative 
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Cluster extension 1 (in units of r^) 
Fig. 2. Influence of the finite extension on the Compton parameter, the surface brightness and the 
Hubble constant assuming f3 = 2/3, r^ — r^ ^ and a spherical cluster. The relative errors 
which are shown are defined in the text. 



errors e^™* on the Compton parameter and e^ "* on the surface brightness are defined by the 



expressions: 



ftm _ y(oo) -yjl) , Um _ Sxioo)- Sxjl) 



(14) 



(15) 



"^ 2;(oo) ^ Sx{oo) 

^From Eqs. (M) and (HO) we can easily estimate these ratios which are 
^,„,_ i3,„ (3/2/3 -1/2, 1/2) j,^,^ B^{3P- 1/2,1/2) 

^y 5(3/2/3-1/2,1/2) '"^"^ ^^ 5(3/3-1/2,1/2) ' 

For a spherical density profile, a line of sight through the cluster center and an infinite extension, 
we obtain from Eq. (p^: 

5x(oo)-^'- '^ 5(3/3 -i,i) ■ ^''^ 

We introduce 6c = Tc/Da the angular core radius, where Da is the angular diameter distance of 
the cluster: 



n C c goz+ go-1 Vl + 2go2:-l 

Da^ -^ where C = — . 

Ho q^ l + z)2 



(17) 



With Eq. (Hq) we can estimate the Hubble constant Ho{oo) for an infinitely extended cluster by: 



Hoi^) = X'T!P 



sx{^), [5(1/3 -ii; 



with A = AC 



(18) 



y2(oo) 5(3/3 -i,i) 

and for a finite extension I, we get instead 

TT (,\ _ yrr.3/2 ^ X jl) ^ [B (f/3 - ^ , ^) ^ B,n (f /3 - ^, ^)] 

mn-X^eo y,^^6c [5(3^_1^1)_5„^(3^_1^1)] • 

Obviously, Sx andi/'^ are observed quantities and thus the ratios 5x(oo)/y^(oo) and S'x (0/2/^(0 
are in the following both set equal to the measured value ( Sx/y"^ ) • This way, for our exam- 

V / obs 



(19) 
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pie, we get for the relative error for the estimation of the Hubble constant, due to assuming an 
infinite extension rather than a finite one, the expression: 



ftm _ Hajoo) -fl"o(0 _ 



B(3/3 



2' 2) 



[BJlP- 



— Bmi^P ■ 



1 1)1 2 
2' 2^J 



Ho{^) ^ B^{^(3 



h 5) [-^(3/3 ^ 5' i) ^ 5m(3/9 - 5, 5)] 



(20) 



In Fig. (Ij) we plot the relative error as a function of the finite extension / for a spherical 
isothermal /3 = 2/3-model. For I -^ 00 the relative errors on the Compton parameter, the X-ray 
surface brightness and the Hubble parameter become, of course, negligible. For instance, for a 
cluster with a finite extension of about 10 times r^ the relative error with respect to the assump- 
tion of an infinite extension is only about 6 % for the Compton parameter. For the X-ray surface 
brightness, the relative error due to the finite extension is much smaller, for instance an error of 
about 4% is obtained if the cluster has an extension of only about 2 times r,,. Nevertheless, for 
clusters with a small extension (3rc) with respect to their core radius, the error in e^™* becomes 
quite substantial (^ 20%) for the Compton parameter. 



These estimations are in accordance with the results of Inagaki et al. ( 1995 ). The net effect when 
one considers infinite clusters is thus to overestimate the value of the temperature decrement and 
the X-ray surface brightness. 

The influence of the finite extension on the Hubble constant given in Eq. (|o|) is larger. Indeed, 
the Hubble parameter is overestimated by almost 20%, when considering for instance a cluster 
with infinite extension as compared to one with an extension of Ivc (cf. Fig. 0). 



3.2. Polytropic index 



Although the isothermal distribution is often a reasonable approximation of the actual observed 



clusters, some clusters do show non-isothermal features. Henriksen & Mushotzky (1985) have 
suggested that the isothermal model cannot be consistently applied to gas distributions in clus- 



ters. Indeed, Markevitch et al. (1998) find that the temperature profiles of some clusters can 



be approximately described by a polytrope. More recently, Ettori et al. ( [2000| ), with a combined 
analysis of the BeppoSAX and ROSAT-PSPC observations, showed that a polytropic profile with 
index 7 = 1.16 ± 0.03 fits the temperature distribution of the cluster A3562 very closely. 
As a consequence, within the virial regions of typical clusters of galaxies, the gravitating mass, 
the gas mass and the gas fraction can vary quite substantially, compared to that obtained assum- 
ing an isothermal profile (Ettori & Fabian 1999| , Ettori 2000 Ettori et al. 2000). It is interesting, 
therefore, to investigate the variations due to a polytropic equation of state on SZ quantities. 
A polytropic thermal profile has the following form: 



17-1 



T = T 



(21) 
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where the subscript o denotes the values in the cluster center and 7 is the polytropic index. The 

isothermal profile is obtained by setting 7 = 1. 

Then the Compton parameter and surface brightness are given by 



y 



poly _ 2 






^\/ -^ eo^t^n^c 



qpoly _ 

^ 27r(l + z)3 Jo \n. 



dr„ 



dry 



(22) 



(23) 



Similarly to the calculations in Section 2, we find for yP°'y and S^ ^: 




y'""Hr^,r.)^ - 


rr /- / 2 2 \ -1/37+^ 
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.(!,,_-)_.„,(!,,_-)] 
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Sfyir,,r^)^ - 


47r(l + z)3 1,^+^2 + ^2; 
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/3 /7 3\ 1 1\ /3 / 
«(2''(i + 2)-2'2)-*»(2''( 


'7 3\ 
,2 + 2J 


---11 
2'2A 



(24) 



(25) 

The relative error by considering a polytropic profile (index poly) as compared to an isothermal 
profile (index iso), can be expressed as follows, assuming for simplicity an infinite extension of 
the cluster, 



^poiy _ y± 



Vpoly 



yiso 


= 1- 


VCi/ VCa 


for the Compton parameter and 


poly WXJiso ~ WXJpoly 


{Sx\so 


= 1- 


VCi/ VCa 



(26) 



1/3(7-1) 



5(1/37- 



1 1^ 

2' 2; 



V 2 / 2 ' 2 ^ 



(27) 



|/3(7~1) 



5(|/^(i + i)-ii) 

i?(3/3-i,i) 



for the surface brightness. 

The corresponding error for the Hubble constant turns out to be 



poly _ i^o)iso ~ (^o)po 



ly 



'H 



1 



[B{lPl~hh)]'B{3P 



1 1^ 

2' 2> 



[i3(i/3-i,i)]^i3(i/?(i + i)-i,i) 



(28) 



(■f^o)jso 

when taking measurements along the line of sight going through the cluster center. The relative 
errors defined above for polytropic indices between 1 and 2 are shown in Fig. (pt). Already, a 
small deviation from the isothermal case (7 = 1) leads to significant relative errors and thus to 
quite different values of the observable quantities. For instance, assuming an isothermal instead 
of a polytropic profile with index 7 = 1.2 leads to an overestimation of 20%, 4% and 34% for 
the quantities y, Sx and Ho, respectively. 
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Fig. 3. Relative error e^"'^ (for ATg^), 6^°'^ (for Sx) and e^°'^ for the Hubble constant between 
a polytropic and an isothermal profile. The line of sight is taken to go through the center of the 
cluster, which is assumed to have a spherical profile with infinite extension (/3 — 2/3-model). 



33. Geometrical effect 

In order to compare the relative errors induced by geometrical effects, i.e. the difference between 
ellipsoidal and spherical geometries, we choose the parameters of the cluster for both geometries 



such that we get the same value for the emission integral EI, defined by Sarazin (1988) as: 
EI == /" nl dV, (29) 

where V is the volume of the cluster We consider, for simplicity, the case where the various 
integrals are approximated by assuming an infinite extension: after some algebra (see Appendix 
A for the change of variables) we get the following expression for the emission integral for a 
spherical geometry: 

r(:\B- I) 

(30) 



,3/2 2 3 
''eo' c 



r(3/3) 

and for an ellipsoidal geometry: 

Thus the equal emission condition Elgph — EI^u requires Ci (,2 Cs = 1- 



(31) 



In the following we will consider the two axisymmetric cases prolate (cigar shaped) with sym- 



metry axis Tx, thus C2 = Cs = vVCi' and oblate (pancake shaped), where the symmetry axis is 
given by r^, and thus C2 = Ci and C3 = 1/Ci ■ 



We introduce ef,*^"™ and e 



^g as the geometrical relative error between ellipsoidal and spher- 



ical geometry for the Compton parameter and the surface brightness, respectively, 

geom _ Vsph ~ Veil , geom _ \'-'X)sph ^ \'^X fell 



Vsph 



yell 1 neoni 

and e^ = 



iS- 



(32) 



X, 
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Fig. 4. Geometrical relative error e^*^"™ for y between an axisymmetric ellipsoidal and a spherical 
geometry (assuming infinite extension, isothermal profile and (3 ~ 2/3) for three different lines 
of sight, parametrized in units of Tc. The two ellipsoidal shapes prolate and oblate are indicated. 



and thus find[| 



geom 



I-C2 



= 1-C2 




l/3+i 



-3/3+ i 



and 



(33) 



(34) 



In Figs. (0) and (pb we have plotted the geometrical relative error as a function of the flattening of 
the profile Ci /Cs for the Compton parameter e^'^""'' and the surface brightness Cg^"™', respectively. 
We consider three different lines of sight, {rx, r^) = (1,0), (1,1) and (0,1), given in units of the 
core radius Vc- 

Fig. (0) shows, that the line of sight (1,0) in the prolate case leads to almost the same error of 
the Compton parameter as the line of sight (0,1) in the oblate case. For a flattening of 50% (i.e. 
C1/C3 = 1.5) we get an overestimation of about 2% by using a spherical cluster shape instead 
of an elUpsoidal one. The other lines of sight lead in the prolate case to overestimations of up to 
19%, while in the oblate case the Compton parameter is underestimated by almost 22%. 
Due to a quadratic dependance on the density profile, the surface brightness is more affected by a 
flattened shape. The Une of sight (1,1) shows in both cases (prolate and oblate) an overestimation 
of about 8%. Towards (0,1) the surface brightness is overestimated in both cases by up to 30%, 



' Remember that the set of coordinates (r^;, rj,, r^), as well as the half axes of the ellipsoid ((^1, (^2, Cs), 
are given in units of the core radius re- 
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Fig. 5. Geometrical relative error gl"^"™ for Sx between an axisymmetric ellipsoidal and a spher- 
ical geometry (assuming infinite extension, isothermal profile and (3 ~ 2/3) for three different 
lines of sight, parametrized in units of Tc- The two ellipsoidal shapes prolate and oblate are 
indicated. 



while the view towards (1,0) results in an underestimation of 24% for a prolate deformation and 
almost 38% in the oblate case. 

In the case of an infinite extension approximation from Eq. ( [l2| ) we get 

|^^ATi'^C,r.x(l + ^ + ^) ^ ,l^^_lf (35) 

Again, the ratio y'^ /Sx ratio is set equal to the observed value and leads then to the relative error 
on the Hubble constant: 

Fig. (^ shows in the axisymmetric elUpsoidal cases, prolate and oblate, the influence of a flat- 
tening of the cluster profile up to 50%. The prolate case results in a systematic overestimation of 
the Hubble constant, by considering a spherical instead of a flattened profile. Depending on the 
line of sight, this error goes up to 22%. Underestimations of up to 33% arise in the oblate case. 
We summarize in Table (jl|) the relative errors on the i/-parameter, the surface brightness and the 
Hubble constant that appear by considering a spherical instead of a flattened ellipsoidal shape 
with, as an example, an axis ratio of 1.2. Negative values indicate that the considered quantity is 
underestimated, whereas positive values that it is overestimated. 

In the general triaxial case (i.e. not necessarily prolate or oblate), estimating the Hubble 
constant is more difficult. We can compute the angular diameter distance Da and thus the Hubble 
constant by measuring the angles 9i and 6*3, defined as the angular ellipsoidal core radii 9i = 



12 
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Fig. 6. Geometrical relative error gS^"™ for H between an axisymmetric ellipsoidal and a spher- 
ical geometry (assuming infinite extension, isothermal profile and f3 — 2/3) for three different 
lines of sight, parametrized in units of re- The two ellipsoidal shapes prolate and oblate are 
indicated. 



Cluster shape 


Line of sight 
(in r-c) 


geom 

(in %) 


geom 

(in %) 


geom 

(in % ) 


prolate 


(1,0) 


0.4 


-11.7 


11.1 




(1,1) 


4.5 


1.7 


7.3 




(0,1) 


8.8 


14.3 


2.9 


oblate 


(1,0) 


-9.4 


-16.1 


-3.2 




(1,1) 


-3.5 


1.7 


-9.1 




(0,1) 


0.4 


12.4 


-13.3 



Table 1. Relative errors on the y-parameter, the surface brightness and the Hubble constant are shown. The 
flattened ellipsoidal is supposed to have an axes ratio of 1.2. Negative numbers indicate underestimations, 
whereas positive ones overestimations. 



C^iTc/ Da and 6*3 = Ca^c/i'A, but we have no observational access to the corresponding angle 
62- On the contrary, in the case of a spherical profile it is sufficient to measure, as we have seen 
in section (3.1), the value of y and Sx towards the center of the cluster in order to evaluate the 
Hubble constant. We do not discuss further the geometrical relative error for a general triaxial 
cluster 
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3.4. Projection effect 



Ruiz ( 1976| ) and Stark ( 1977 ) have discussed the projection onto the sky of luminosity distribu- 



tions which have an ellipsoidal form. This problem has also been treated by Fabricant, Rybicki 



& Gorenstein ( 1984 ). The projection effect is expected to broaden the measurements of the tem- 



perature decrement and the surface brightness as mentioned by Cooray ( 1998 1. 
In order to evaluate the projection effect we assume an ellipsoidal geometry. As an example, we 
take the case of an infinite cluster extension with isothermal profile and a prolate form. 
The profile of the cluster is supposed to be rotated by an angle 6 around the r^ axis. For compar- 
ison, we consider a prolate profile with a major half axis along r^, which has the same projected 
image on the sky, i.e. on the (r^, r2)-plane. We then compute the difference on the resulting y- 
parameter and surface brightness between these two profiles. The structure integrals IP^°J and 
jrot ^g given in Appendix B, where we have computed the y parameter and the surface bright- 
ness of a cluster both for a rotated coordinate system, with an angle 9, and by considering its 
projection onto the sky. 
The relative error can be quantified as 

Tproj Trot 

(37) 



/ jproj 

It is a pure geometrical effect and thus the same for the y parameter, the surface brightness and 
the Hubble constant. 

In Figs. (0) and (H) we have plotted the relative error due to the projection effects as a function 
of the rotation angle 6 and the axes ratio Ci/Ca of the ellipsoid. 

The maximal 9 is assumed to be 7r/4 ^ 0.785, for which we find an underestimation of almost 
17%. The influence of the axes ratio turns out to be less than 5%. 

4. Discussion and conclusions 

We have calculated the relative error caused by assumptions regarding finite extension, a poly- 
tropic temperature profile, ellipsoidal geometry and projection effects, on the measurements of 
the X-ray surface brightness, the SZ temperature decrement and the determination of the Hub- 
ble constant. Although the X-ray data have improved dramaticafly in the last decade, it is still 
difficult to determine the internal structure of clusters from X-ray imaging alone, because such 
images supply only projected temperature and surface brightness information, without further 
indications of the internal gas dynamics. Nevertheless, recent observations show indirectly that 



many clusters are still dynamically evolving (see Mohr et al. 1995). 



Cooray (2000) has discussed intrinsic cluster shape, in particular considering axisymmetric mod- 



els such as oblate and prolate elUpsoids, using the Mohr et al. (y_995|) cluster sample. Their study 



shows that clusters do indeed have aspherical profiles, which are more likely described as prolate 
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Fig. 7. The relative error, when comparing a projected and a rotated prolate-shaped cluster, is 
shown as a function of 9. The axis ratio is fixed to be Ci/Ca — 1-5- 
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Fig. 8. The relative error between a projected and a rotated prolate-shaped cluster profile is shown 
as a function of the axes ratio Ci/Cs- The rotation angle 6 is fixed at n/A. 



rather than oblate ellipsoids. Nevertheless, Mohr et al. (1995) remarked that they cannot rule out 
the possibility that clusters are intrinsically triaxial. 

Pierre et al. (199q) studied with ROSAT the rich lensing cluster Abell 2390 and determined 
its gas and matter content. They found that on large scales the X-ray distribution has an elliptical 
shape with an axis ratio of minor to major half axis of Ci /Cs ~ 1 .33. Using our results we see that 
this corresponds to a relative error in the y parameter of up to 10%, depending on the line of sight 
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and the shape of the cluster (prolate or oblate, see Fig. m. The surface brightness measurements 
lead to errors of up to 25% (see Fig. Bb and thus the Hubble constant is overestimated by about 
23% (see Fig. |). 

An unresolved temperature gradient in the gas affects the gas profile and thus the total mass 
derived assuming hydrostatic equilibrium. If such a gradient is present, the true temperature in 
the central region may be higher than the emission-weighted temperature generally used. As an 



example, Grego et al. ( ZOOO ) observed in Abell 370 a slow decline of the temperature with radius. 
The temperature falls to half its central value within 6-10 core radii. This temperature profile can 
be approximately described by a gas with a polytropic index of 7 — 1.2, which in itself is already 
an important modification with respect to an isothermal profile and could lead to a relative error 
e^ ^ of 37% in the evaluation of the Hubble constant (see Fig. Bj). 

Furthermore, the optical and X-ray observations of this cluster show a possible bimodal mass 
distribution. Thus, the combined temperature and geometry effects must be taken into account to 
obtain reliable values for such parameters as the gas and total matter content. A similar polytropic 



index (7 = 1.16) has also been found for Abell 3562 (Ettori et al. 200C ). 

Cooling flows in galaxy clusters can substantially change the temperature profiles, especially 



in the inner regions. Schlickeiser ( 1991 ) and Majumdar & Nath (2000 1 have investigated the 
changes induced by a cooling flow in the temperature and density profiles, and their implications 
on the SZ effect. We notice that for a polytropic distribution, the density profile can still be 
well approximated hy a (3 profile, whereas for cooling flow solutions the density becomes quite 



different. For example, Vikhlinin et al. (1999) showed that outside the cooling flow region, the 
/3-model describes the observed surface brightness closely, but not precisely. In this context. 



Majumdar & Nath ( 2000 1 found that the presence of a cooling flow in a cluster can lead to an 
overestimation of the Hubble constant determined from the SZ decrement. 



Recently, Mauskopf, Ade, Allen et al. (2000) determined the Hubble constant from X-ray 
measurements obtained of the cluster Abell 1835 with ROSAT and from the corresponding mil- 
limetric observations of the SZ effect with the Sunyaev-Zel'dovich Infrared Experiment (Suzie) 
multifrequency array receiver Assuming an infinitely extended, spherical gas distribution with 
an isothermal equation of state, characterized by /3 = 0.58 ± 0.02, T^o — O.Sl^'g keV and 
Ueo = 5.64ti[o2 X 10~^ cm^^, they found a value of Ho ^ ^^tH km s^^ Mpc^^ for the Hub- 
ble constant. In Fig. (m we show the influence of geometry and of assumptions of finite extension 
on this result using the same input parameters. Fig. (9a) shows that for a spherical geometry. Ho 
displays a strong dependence on the cluster extension. Fig. (9b) gives the value of Ho assuming 
an infinite extended ellipsoid shaped cluster (instead of a spherical geometry), as a function of 
its axis ratio Ci /Ca- 
in summary, we see that it is crucial to know the shape of a cluster and its temperature profile. 
For this problem, the new X-ray satellites have the necessary spatial and spectral resolution to 
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Fig. 9. The Hubble constant derived from the data of Mauskopf et al. (2000). Fig. (9a) shows the 
influence of finite extension, while Fig. (9b) gives the value of Ho assuming an axisymmetric 
ellipsoidal geometry. In the latter case, oblate or prolate geometry give the same value of Ho 
when taking a line of sight through the cluster center, as is assumed here. 



remove the effects of contaminating sources in the field and to measure the spatial variation of the 
cluster temperature. In this context it will be better for future studies to focus on nearby cluster 
samples, which are less subject to observational selection effects (as mentioned by Roettiger et 



al., (1997)) 
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Appendix A: 



Structure integral for the y parameter 

The structure integral for the SZ effect is given by the expression 

Isz{rx,rz)^2 —dry, (Al) 

Jo '^eo 

where Ue is the density profile of the electrons, Ueo the central density of the cluster and I the 
maximum extension of the hot gas on the line of sight ry, in units of the core radius re- With an 
ellipsoidal density profile, we obtain 



Isz{rx-,r^) = 2 



^A <T»^ „2 

'- ' /■2 ^ /■2 ^ /■2 
Si 1,2 ^3 



\0 



dr„ 



,y, (A2) 

where r^ , Q and I are given in units of the core radius Tc- With the definition of a function G and 
transforming the variable of integration from Vy to ^ such that 

2 2 2 

G(o = (1 + 0^'^'' with ^ = !;i + ^ + !;i , 

'.I S,2 S3 

we obtain after some algebra the structure integral 

f-'?+(VC2)' „2 



G{£.) {£, - fl)-^'^ di , where r; =. -| + ^ . 

Si 



r 

z 
S3 



Finally, with a last change of variable. 



we get 
with 



1 



77 



1+e ' 



1 



n 



i+^+a/C2)2 

The structure integral turns out to be 

where we introduced the Beta and the incomplete Beta functions defined by 

r(o)r(i.) 



^' 1^-5-5) -«™(l'^-?5 



and 



B{a,b) 



Bmia^b) 



T{a + b) 
Tm{a)Tmib) 



Tmia + b) 
with the Gamma and the incomplete Gamma-functions F and F^, respectively. 

X-ray surface brightness 

The structure integral for the X-ray surface brightness is given by 



^Sx 



I r 



dr„ 



With the same transformations, as given above, we get 

IsAr.,r,) = C2 (1 + vr'^+i B Up -\^^- ^n 



^ 2'2 



(A3) 



(A4) 



(A5) 



(A6) 



(A7) 
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Projection effects on the y parameter 

The rotation in the {tx, rj,)-plane around r^ with an angle 9 leads to the density profile: 



'^e vxiTyiTzi 



1 



[vx cos9 + Ty sin6')^ {vy cos9 — r^ sin6 



CI 



'=.3 



-3/3/2 



(Bl) 



Cl S2 

The rotation angle 9 is the angle between the major half axis of the rotated ellipse in the [r^ , ry ) 

plane and the r^j-axis. The line of sight is taken to be along the Ty-axis. 

To investigate the projection effects we will compare a rotated cluster with respect to the r^ axis 

with its projection on the irx,rz) plane, corresponding to the sky plane. We assume an infinite 

extension and an isothermal profile. 

The rotated structure integral for the y parameter turns out to be 



im^ 



1 



[tx cos9 + ry sin6 



(ry cos9 — rx sin6 

(I 



S3 



-3/J/2 



dry .(B2) 



Cl S2 

After some algebra and with the same kind of variable changes as in Appendix A, we get (as- 



suming an infinite cluster extension for the structure integral): 



with 



1 



s2 






1 1 

2' 2 



(B3) 



m 



C?cos26' + Clsin^9 Q 



On the other hand, the projection of this cluster on the observed sky plane in the infinite cluster 
extension case leads to the density profile 



1 



-3/3/2 

'JL 

C( Cl 

where (i is the maximum value that we get along the rx axis in units of re 



^^"^rx.ry.rz) = rie 



'»3 



(B4) 



Ci^^Cico^^9 + Cisin'9, 

moreover, 1^3/(2 = C3/C2 = 1 (prolate). 

The projected structure integral for the y parameter turns out to be 

with 



^sz\rx,r.,9) ^ C,2{l + rie)-^-"+-- B [-P - -,- 



3 



1 1 



(B5) 



c? 



C3^' 



X-ray surface brightness 



For the structure integral of the X-ray surface brightness we find in the case of the rotated cluster 
rSi{r.,,r,,e)= , ^: __ (l + „)-w<B|3,3-^.^| , (B6) 






and for the projected case 



^.■.-^f^HB[iP-\,\ 



(B7) 
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